Excitations of the field-induced soliton lattice in CuGeOs 
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Here we report the first inelastic neutron scattering study of the magnetic excitations in tfie in- 
commensurate high-field phase of a spin-Peierls material. The results on CuGeOs provide direct 
evidence for a finite excitation gap, two sharp magnetic excitation branches and a very low-lying 
excitation which is identified as a phason mode, the Goldstone mode of the incommensurate soliton 
lattice. 



A one-dimensional (ID) spin ^ Heisenberg antiferro- 
magnet is unstable with respect to dimerization as well 
as long-range antiferromagnetic order. Coupled to a 
three-dimensional (3D) phonon field, it can undergo a 
second-order phase transition at a finite temperature 
into a dimcrized (spin-Peierls) ground state with total 
spin Stot — '^iSi — (D-phase). The dimerization 
is supported by antiferromagnetic next-nearest-neighbor 
exchange in the ID chain direction. Antiferromagnetic 
(AF) interchain couplings prefer a Neel-type ground state 
and compete with dimerization. 

In the D-phase of a spin-Peierls material, the lowest 
magnetic excitation is an isolated triplet branch {Stot = 
1). This excitation is a bound pair of domain walls with 
respect to the dimer order parameter. The two domain 
walls, or solitons, are created by breaking and delocaliz- 
ing one dimer bond. The binding energy depends on the 
magnetic and elastic interchain interactions. The energy 
of the triplet is finite over the entire Brillouin zone and 
achieves its minima at wave vectors k = and k = - 

c 

(the AF zone center, c being the average distance be- 
tween two spins in the chain direction.) In a magnetic 
field, the Stot = 1 branch becomes Zeeman split, and a 
magnetized incommensurate (IC) phase is entered above 
a critical field He, which corresponds approximately to 
the field where the lowest mode should soften to zero en- 
ergy. In the IC-phase, the magnetization is generated by 
introducing more and more domain walls into the dimcr- 
ized ground state. The unpaired spin at the center of 
a domain wall aligns parallel to the magnetic field (as- 
sumed to be II to z). The transverse parts of the spin 



citation gaps, because at each given magnetic field the 
lattice, and hence the intrachain exchange, adapts to the 
magnetization and the number of solitons. Two excita- 
tion branches with polarization perpendicular to H are 
predicted, A± , corresponding to an increase or decrease 
of the total spin by 1, with minima at A: = and 
k ~ where m is the magnetization per spin in units of 
/is [|^ . A longitudinal branch Aq should also be present 
with minima at fc = and fc = — m-), and with 

A™™ > i(A!p" + A!!^™) |]. At^'fc = f this branch 
should have an energy qhbH 0. The precise value and 
field dependences of the modes A± and Ag depend on 
the details of the assumptions used in the models [^-|| . 

CuGeOs is the first spin-Peierls compound where large 
single crystals are available thus allowing the com- 
plete phonon and magnetic excitation spectra to be mea- 
sured using neutron scattering techniques. The IC phase 
is also accessible, since magnetic fields up to 14.5 T have 
become available for neutron experiments, well in excess 
of the critical field in CuGeOs of hqHc ~ 12.5 T. Our 
recent neutron scattering study of the magnetic soliton 
lattice in the IC-phase reveals a static magnetic mod- 
ulation with the same period as the distortive modula- 
tion, and in particular the superlattice Bragg peaks at 
(^,1, \±5ksp) were found to be almost entirely magnetic 
in origin. Outside the critical field region, the magnetic 
soliton width, and the amplitudes of the staggered and 
uniform magnetic modulation are in reasonable agree- 
ment with the predictions of field theory. 

In this letter we report on the magnetic excitation 
spectrum in the IC phase, measured both close to the 



Hamiltonian, SfS^ , delocalize the solitons. This results superlattice peak (5,1,5 ± Sksp), and also at the AF 



in an equal spacing of solitons, and thus leads to an IC su- 
perlattice of distortive and magnetic solitons. Their finite 
width implies a staggered magnetic polarization close to 
the soliton center. 

For a spin-Peierls system, the magnetic excitation 
spectrum in the IC-phase is expected to have finite ex- 



zone center (0, 5, 1), where the dispersion of the magnetic 
excitations in the D-phase has its minimum. The exper- 
iments were performed at the cold triple axis spectrome- 
ter FLEX at BENSC, Berlin, using the vertical cryomag- 
net VM-1. Single crystals of CuGeOa of 0.34 cm^ and 
0.49 cm-^ were aligned in the (0, fc,^) and {h,2h,£) scat- 
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tering planes, respectively. The collimation was guide- 
60'-60'-60' with either fc/=1.3 A fixed and a nitrogen- 
cooled Be-filter in fc^, or fc/ = 1.54 A and the Be-filter 
in kf. The critical fields at (0,1,^) and (^,1,^) dif- 
fer (12.58 T and 12.31 T) because of the different field 
directions and hence g-factors (2.152 and 2.199) in the 
respective scattering geometry. 
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FIG. 1. Field dependence of the triplet excitation in 
CuGeOa at (0,1, i), 2 K. Solid lines: triplet splitting 
ga^J,BHStat, 9a = 2.152. The critical field is fioHc = 12.58 T. 
The inset shows the high-field region in more detail. 

Figure |l| shows the magnetic field dependence of the 
excitations at the AF zone center (0, 1, ^). With increas- 
ing field, the lowest triplet branch decreases to about 
0.4 meV but never becomes completely soft. At the crit- 
ical field, two new sharp excitations appear at ~1 meV 
and ~0.5 meV. The 1 meV mode shows barely any field 
dependence between 12.5 T and 14.5 T, while the energy 
of the lower mode increases towards the higher mode. 
The excitations at (^, 1, ^), close to the magnetic Bragg 
peaks at (i, 1, ^ ± Skgp), display the same field depen- 
dence as those at (0, 1, i), apart from an overall shift 
towards higher energies (1.2 meV and 1.6 meV) which 
is entirely accounted for by the ferromagnetic interchain 
interaction along a*. 

In Fig. ^ we show high-resolution scans of the low- 
energy region at a wave vector of (i,l,i) for selected 
magnetic fields. At 12 T, the widths of the middle and 
the lowest triplet excitation are resolution limited. Both 
broaden significantly at the critical field (about 12.3 T) 
and then disappear. In a small field region, the new 
sharp IC-phase excitations at 1.2 meV and 1.6 meV co- 
exist with the broadened D-phase modes (see scan at 
12.5 T), confirming the first-order character of the tran- 
sition which has also been observed with other methods. 
No evidence was found for a third excitation at an energy 
of gfisH at any of the fields investigated. 

The spectral weight of the two new modes in the IC- 
phase is concentrated in a very small region extending 
from the commensurate position £ = ^ to the incom- 



mensurate one at £ = ^ ± Skgp. The upper part of Fig. 
^ shows the dispersion of the lowest two excitations at 
(i,l,i) along the chain direction c*, at T = 2 K and 
B = 14.5 T, where 6ksp — 0.015. The dispersion (solid 
lines) appears significantly flatter than the dispersion cal- 
culated from the zero field spin-wave velocity (dashed 
lines). Neither at (0, 1, ^) nor at (^, 1, ^) do the disper- 
sions of the two sharp modes exhibit minima at the IC 
positions £ — ^ zt 6ksp- 
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FIG. 2. Energy scans at (i, 1, |), T=2 K. Solid lines: fits 
to Lorentzian-type profiles including the Bose factor. The 
critical field is ^qHc = 12.31 T. 

A search was also made for excitations emanating from 
the magnetic IC Bragg peaks at (5, 1, 5 ± Sksp). This 
revealed a mode at -0.29 mcV (Fig. |), at T = 2 K and 
B = 14.5 T. The energy of this mode (at the respective 
IC wave vector) changes only little between 13 T {Skgp ~ 
0.011) and 14.5 T {Sk^p = 0.015). The dispersion of this 
mode was studied at B = 14.5 T and is plotted in the 
bottom part of Fig. ^. Its spectral weight is concentrated 
at the IC wave vector (i, 1, i ± Skgp), the minimum of 
its dispersion. 

In interpreting the data the first thing to establish is 
whether the three modes shown in Fig. |^ are magnetic 
or nuclear in origin. It is known from our earlier study 
that the satellites at (5,1,^ ± Sksp) are predominantly 
magnetic (by a factor of ^ 40) m, and it follows that 
the low-lying mode emanating from the satellites must 
also have a predominantly magnetic character. For the 
higher-lying modes it was observed that their intensity 
decreased at equivalent k but larger total wave vector 
transfer. They too can therefore be ascribed to be mainly 
magnetic in origin. It should be noted, however, that a 
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structural contribution cannot be excluded, and indeed 
a mixing of magnetic and elastic degrees of freedom is 
expected in a spin-Peierls material. 
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FIG. 3. (i, 1, ^)-dispersion of the lowest magnetic ex- 
citations in CuGeOa at 14.5 T, 2 K. Solid: disper- 
sion fit to V ^min + (^-^min)^, where £min = 5 fot 

A± and ^ ± Sksp for the 0.29 meV excitation (phason), 
and AA^,A_,ph = 41(5), 48(5), 32(2) meV/r.l.u.. The mag- 
netic Bragg peaks are at i ± Sk^p with Sksp — 0.015. 
Dashed: dispersion assuming the H — spin-wave velocity, 
A = 96.3meV/r.l.u.. 
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FIG. 4. Field dependence of the minimum energy of the 
phason, energy scans at 12 T, 13 T and 14.5 T at the respec- 
tive modulation vector (i, 1, i -I- Sksp). 

Calculations of the excitations in the IC-phase are 
reported for three cases, with dimerization parameters 
e = tJ^t^ = 0.014 §, e = 0.14 |], and e = 0.4 
[||. Poilblanc et al. and Schonfeld et al. @ also in- 
cluded frustrating ncxt-ncarest-neighbor intrachain ex- 
change a — -rr-f — — 0.24 and a = 0.36, while in the cal- 



culations of Yu et al. |2| only nearest-neighbor intrachain 
exchange was considered. Irrespective of this large spread 
of the parameters e and a the spectral weight is concen- 
trated at the commensurate position ^ = \ @J§li where 
A± have minimum excitation energy [Q. This agrees 
with the observed upper two modes at 1.3 meV and 
1.7 meV (14.5 T, Fig. 3), which we therefore identify with 
the modes A± . The longitudinal mode Aq should have its 
minimum at qc = \ ±6ksp with A™" > i(A™™ + A™™) 
11 and g^iBH at gc = ^ [§. Therefore the 0.29 meV 
excitation cannot be identified with Ap. A weak longi- 
tudinal contribution at higher energies could be hidden 
by the transverse excitations. However, the weakness or 
absence of the longitudinal mode is probably related to 
the significant interchain interaction in CuGeOa. Cal- 
culations which imply significant spectral weight of 
the longitudinal mode do not include interchain inter- 
action, and it is known that longitudinal modes usually 
become stronger with increasing one dimensionality of 
the exchange [p|jlC|]. 

While the field dependences of the modes A± could 
in principle be used to place further constraints on the 
values of the coupling parameters, this turns out to be 
difficult in practice. This is mainly due to two reasons, 
to the fact that the values of e used in the calculations 
lead to gaps much larger than reported here, and to 
the fact that the magnetization range studied starts at 
larger values than the maximum magnetization m=0.015 
/is in our experiment. It is, however, perhaps worth 
noting that the calculation with a finite and large value 
of a Q predicts that A+ should increase with increasing 
field, while A_ should be field independent, in qualitative 
agreement with the data shown in Fig. |l|. In this sense 
our data are consistent with the large value of a needed 
to explain various experimental facts |lll|-p^, and with 
the notion that the spin-Peierls transition in CuGeOa is 
driven by hard phonons [jisl-^. 

We now address the origin of the observed excitation 
at 0.29 meV. Its field dependence is too small for an 
optic mode arising from a folding back of the Brillouin 
zone. Such an optic mode should shift with 5ksp by about 
0.1 meV from 13Ttol4.5T.Afittoa Lorentzian hne- 
shape including the Bose factor for T = 2 K yields the 
finite energies shown in Fig. ^. The same gap energy 
(2.3±0.2 cm~^ « 0.29 meV) was observed in a Raman 
experiment |]l8| , which mainly probes the structural ex- 
citations, pointing to a magneto-structural origin for this 
excitation. We attribute the mode at 0.29 meV to phase 
oscillations of the soliton lattice, so-called phasons. The 
IC long-range order in the high-field phase breaks the 
quasi-continuous translation symmetry of the soliton lat- 
tice with respect to simultaneous shifts of the soliton 
centers along the chain axis. The Goldstone modes cor- 
responding to this symmetry breaking are phase oscil- 
lations of the magnetic and distortive soliton structure. 
The corresponding eigenvectors carry magnetic and elas- 
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tic amplitudes. Their dispersion close to the IC superlat- 
tice Bragg peak is expected to be given by the spin- wave 
velocities of the dimerized phase, except close to the crit- 
ical field ||f^ . Indeed the experimentally observed disper- 
sion along £ is of the same order of magnitude as that of 
the transverse magnetic excitations (Fig. and is flatter 
by a factor of 3 than the zero-field spin-wave velocity. 

If the mode at 0.29 meV is indeed a phason, it re- 
mains to be explained why it has a gap. For a spin- 
Peierls system with soft phonons, a first-order transi- 
tion into the high-field phase could proceed via a series 
of lock-in transitions, which implies a finite phason gap 
|^,|l|. No soft phonons have been found in CuGeOs, 
and a number of other possible mechanisms need to be 
explored, including: impurities, discreteness of the lat- 
tice, and k dependent spin-phonon coupling. A pin- 
ning to impurities would create a disordered soliton lat- 
tice, in contradiction to our earlier structural study 0. 
The discreteness of the lattice leads to a phason pinning 
energy of the order 0.003 meV which is much too 
small. However, this estimate neglects interchain cou- 
plings, and does not take into account the non-sinusoidal 
magnetic modulation |Q, both of which may enhance 
the phason gap in CuGeOs. In the absence of soft- 
phonons, lock-in phenomena in CuGeOs could be pro- 
duced by a A: dependent spin-phonon coupling. Schonfeld 
et al. Q introduce a fc-dependent effective elastic con- 
stant K{k) — X^A ^-^(^)/5a(^)' where VL\{k) are the 
dimerization phonons and gx(k) the respective coupling 
constants. The relevant phonons in CuGeOa do not dis- 
play a softening in the dispersion at the commensurate 
wave vector kgp p^ , but the spin-phonon coupling in- 
creases towards fcgp {K{ksp) ~ 4:K{0) |2^). This prefer- 
ence for a commensurate structure enhances the influence 
of the discreteness of the lattice, and may lead to a finite 
phason energy and a commensurate wave vector, instead 
of gapless phasons and a truely IC modulation. 

The hysteresis observed in ESR experiments in the 
whole IC-phase supports the picture of a lock-in 
of the soliton centers to the discrete lattice. Neverthe- 
less, in the past, several experimental results were taken 
as evidence for gapless phasons: i) an increase of the 
T^-contribution of the specific heat capacity in the IC- 
phase ; ii) the absence of discrete lines in the NMR- 
lineshape; iii) the reduced ratio TOs/to„ of the staggered 
and uniform magnetic modulation amplitude (rris, "t-u) 
derived from NMR experiments. At closer inspection i) 
could be explained by an altered elastic constant in the 
IC-phase rather than by gapless phasons. Phasons with 
a dispersion similar to magnons yield -contributions 
with a' = 1 or 2 corresponding to their ID or at least 2D 
density of states, ii) The NMR lineshape could equally 
well be explained by an unresolved assembly of the ex- 
pected ^^150 (discrete) lines within 0.5 T and hence is 
still consistent with the presence of lock-in phenomena, 
iii) Zero-point fluctuations of gapless phasons should lead 



to the same reduced ms/m„ in a neutron experiment 
where no such reduction was observed. The neutron re- 
sult agrees within 20% with field theory where phasons 
are neglected. This shows that zero-point phason oscilla- 
tions are suppressed and hence indicates independently 
a finite energy gap of the phasons. 

In summary, we have observed the predicted transverse 
magnetic excitations A± in the IC phase of a spin-Peierls 
material which correspond to the creation or destruction 
of a domain-wall pair. The longitudinal magnetic exci- 
tation Aq was not found, probably due to the 3D inter- 
actions. The phason mode of the IC modulated struc- 
ture has also been observed, and is gapped with a maxi- 
mum spectral weight at the position of the IC satellites, 
(^, f , 5 ± Skgp). We argue that the finite energy of the 
phason mode is consistent with all of the experimental 
data available to date. We hope that our data will stim- 
ulate and guide the development of a microscopic model 
of CuGeOs which will ultimately provide a consistent de- 
scription over the entire temperature-field phase diagram 
of this unconventional spin-Peierls system. 
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for funding. 



[2] 
[3] 

[4] 

[5] 

[6] 

[7] 
[8] 
[9] 



[10] 

[11] 

[12] 
[13] 
[14] 

[15] 
[16] 



G. Uhrig, F. Schonfeld, and J. Boucher, Europhys. Lett. 
41, 431 (1998). 



W. Yu and S. Haas, |gond-mat/9909093| (1999). 

G. Uhrig, F. Schonfeld, J.-P. Boucher, and M. Horvatic, 
Phys. Rev. B 60, 9468 (1999). 

F. Schonfeld, G. Bouzerar, G. Uhrig, and E. MuUer- 
Hartmann, Eur. Phys. J. B 5, 521 (1998). 

P. H. M. van Loosdrecht et al, Phys. Rev. B 54, R3730 

(1996) . 

M. Ease, I. Terasaki, and K. Uchinokura, Phys. Rev. 
Lett. 70, 3651 (1993). 

H. R0nnow et al, Phys. Rev. Lett. 84, 4469 (2000). 
D. Poilblanc et al, Phys. Rev. B 55, R11941 (1997). 

M. Steiner and A. R. Bishop, in Solitons, edited by 
S. E. Trullinger, V. E. Zakharov, and V. L. Pokrovsky 
(Elsevier Science Publishers B. V., Amsterdam, 1986), 
Chap. 16, p. 783. 

I. Affleck and G. F. Wellman, Phys. Rev. B 46, 8934 
(1992). 

G. Castilla, S. Chakravarty, and V. Emery, Phys. Rev. 
Lett. 75, 1823 (1995). 

P. van Loosdrecht et al, Phys. Rev. Lett. 78, 487 (1997). 
M. Arai et al, Phys. Rev. Lett. 77, 3649 (1996). 

H. Yokoyama and Y. Saiga, J. Phys. Soc. Jpn. 66, 3617 

(1997) . 

C. Gros and R. Werner, Phys. Rev. B 58, R14677 (1998). 
K. Kuboki and H. Fukuyama, J. Phys. Soc. Jpn. 56, 3126 
(1987). 



4 



[17] G. Uhrig, Phys. Rev. B 57, R14004 (1998). 

[18] I. Loa, S. Groncmeyer, C. Thomsen, and R. K. Krcmer, 

Sol. State Conimun. Ill, 181 (1999). 
[19] S. Bhattacharjee, T. Nattermann, and C. Ronnewinkel, 

Phys. Rev. B 58, 2658 (1998). 
[20] R. de Lima and C. Tsallis, Phys. Rev. B 27, 6896 (1983). 
[21] A. Bruce and R. Cowley, J. Phys. C 11, 3609 (1978). 
[22] M. Braden et al, Phys. Rev. Lett. 80, 3634 (1998). 
[23] R. Werner, C. Gros, and M. Braden, Phys. Rev. B 59, 

14356 (1999). 

[24] W. Palme et al, Phys. Rev. Lett. 76, 4817 (1996). 
[25] T. Lorenz et al, Phys. Rev. B 54, R15610 (1996). 



5 



